Dipolar Drag in Bilayer Harmonically Trapped Gases 
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We consider two separated pancake-shaped trapped gases interacting with a dipolar (either mag- 
netic or electric) force. We study how the center of mass motion propagates from one cloud to the 
other as a consequence of the long-range nature of the interaction. The corresponding dynamics is 
fixed by the frequency difference between the in-phase and the out-of-phase center of mass modes 
of the two clouds, whose dependence on the dipolar interaction strength and the cloud separation 
is explicitly investigated. We discuss Fermi gases in the degenerate as well as in the classical limit 
and comment on the case of Bose-Einsten condensed gases. 

PACS numbers: 03.75.Ss,34.20.Gj,67.85.De,71.45.-d 
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I. INTRODUCTION 



In recent years atomic and molecular dipolar gases 
have attracted a lot of interest since the long-range and 
the anisotropic nature of the interaction is expected to 
give rise to new important features at both the micro- 
scopic and macroscopic level. These include, among oth- 
ers, novel effects in the mechanism of the expansion and 
of the collective modes [l[ and in the structure of the 
superfluid [|| and normal @ phase, new exotic phases 
of crystalline nature (see e.g., the recent work and 
new schemes for quantum computation in the presence 
of optical lattices (5J. Some of these effects, in particular 
those concerning the expansion and the collective modes, 
have been already experimentally observed in magnetic 
dipolar atomic gases |6j. The recent progress in the re- 
alization of gases of electric polar molecules, where the 
effect of the dipolar force is particularly strong, is ex- 
pected to open new challenging frontiers in this area of 
research (see 0-Q and references therein). 



The aim of the present work is to propose a drag ex- 
periment induced by the long range nature of the dipolar 
interaction. We consider an atomic or molecular gas har- 
monically trapped in a double well configuration such 
that the overlap between the two clouds and the corre- 
sponding tunneling effect can be neglected (see Fig. [TJ. 
The only force acting between the two gases is of long 
range nature (here and in the following we assume that 
dipoles are oriented in the direction orthogonal to the 
discs, i.e along the z-th axis of Figure 1) and we study 
how the out-of-phase transverse dipole mode is affected 
by the long-range interaction. Displacing one of the two 
clouds out of its equilibrium position and releasing it, will 
excite both the in-phase (center of mass) and the out-of- 
phase dipole modes. On a time scale fixed by the inverse 
of the frequency difference between the two modes, the 
center of mass motion of the first cloud will be transferred 
to the second one. We call this effect "dipolar dra g" i n 
analogy to the well known Coulomb drag (see e. g.. flop 
exhibited by electrons in uniform bilayer systems 




FIG. 1: Scheme of two not overlapping pancake shaped clouds 
of a dipolar gas. The distance between the centers of mass of 
the clouds is 2L. The clouds are harmonically confined in the 
transverse directions x, y. 



II. DIPOLAR DRAG OF THE 
CENTER-OF-MASS MOTION 

We consider a gas confined by a cylindrically harmonic 
potential: 



v tra P ( x , V, z ) = 7; muj 'i[x 2 + y 2 + X 2 (z ± Zq) Z 
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where 2z is the distance between the minima of the po- 
tential along z, A = uj z /uj± is the ratio between the trans- 
verse and longitudinal trapping frequencies and we con- 
sider pancake configurations, i.e., A > 1. Let x^ being 
the center of mass coordinate along x of the i-th cloud. 
The equations of motion can be written as 



dx\ 
~~dt 

dx2 



u)j_xi + a{x x - x 2 ), 
ljj_x 2 - a{x x - x 2 ), 
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FIG. 2: The motion of the center of mass of the two clouds 
for uj ou t = I.Iojx, with u}±/2n = 200 Hz which corresponds 
to a beating time i — n/(u) ut — = 0.025s. Initially, at 
t — 0, only the cloud 1 is displaced from its central equilibrium 
position. 



where a is the coupling between the two bare center-of- 
mass modes. The eigenfrequencies of the previous equa- 
tion are simply Uj n — uj± , for the in-phase sloshing mode 
and ujout = Wj_a/T + 2a/u> 2 _ for the out-of-phase slosh- 
ing mode. Thus in order to determine a we just need 
to determine the splitting uj out — uj± for the dipolar cou- 
pled system. Once the frequency u out is known, we can 
determine quantitatively the evolution of the system as 
described by Eq. ([3]). In Fig. [2]the motion of the coupled 
clouds for a value of uj out — l.lui± (see Sec. Mil below) 
is shown. The beating of the motion is a direct measure- 
ment of the out-of-phase mode frequency, since the time 
at which the initially displaced cloud stops in the center 
is simply i = Tr/(uj out -tJj_). 

In the following we calculate the frequency uj ut as a 
function of the dipolar interaction strength and the dis- 
tance between the two clouds. We will also discuss how 
the equation of state of the gas affects such a frequency. 
The frequency Lo out was recently calculated by Huang and 
Wu [lH in the case of a magnetic dipolar Bose gas, us- 
ing a technique very similar to the one employed in the 
present work. For this reason we mainly focus on the 
case of Fermi gas. Moreover the Fermi statistics allows 
for an easier realization of cold gases of hethero-nuclear 
molecules carrying an electric dipole moment (e.g., the 
recent experiment [HJ]) so that the strength of the dipo- 
lar force can be much larger. 



The interaction potential between two dipoles (di = d 2 
d) has the standard form: 



V D {f u r 2 ,9) = 



d 2 (l-3cos 2 9) 



\r\ - r 2 \ 



(4) 



where 9 is the angle between d, i.e., the z direction, and 

n - f 2 . 

We study the system by means of the following energy 
functional of the cloud densities ra*, £ = 1, 2 



i=l,2 



(5) 



where we introduce the kinetic energy E l kin = 
ft 2 /(20m7r 2 ) J dr\d>TT 2 ni{r)] b / 3 calculated within local 
density approximation, the potential energy E l trap = 
J df\ni(r)Vt ra ^\ and the intra- and inter-cloud dipolar en- 
ergies 



ni(r)rij(r')VD{r — r')drdr' . 



(6) 



In the energy functional Eq. ([S]) we have not included the 
intra-cloud exchange energy. We safely neglect it, since 
for the pancake-like configurations, which we are mainly 
interested in, the direct term is the dominant effect (see. 
e.g., 0). 

In order to study the center-of-mass oscillations in the 
transverse direction (see Fig. [T]) we consider a scaling 
transformation for the densities of the type rii(x, y, z) — >■ 
7ii{x+ei 1 y, z). From the variation of the energy functional 
© we get as expected two modes. The in-phase mode 
is not affected by the dipolar interaction and has a fre- 
quency lo — uj±_ equal to the harmonic trapping one. Con- 
versely the out-of-phase mode is affected by the dipole 
interaction and is characterized by the frequency 



TV 



dridr 2 V D (fi - r 2 ) 



dni(fi) <9n 2 (r 2 ) 



dx 



dx 



where N is the total atom/molecule number. 

Assuming the two clouds are identical we can get the 
densities of the clouds by means of an easy variational 
gaussian ansatz 



III. TWO COUPLED FERMI GASES WITH 
DIPOLAR INTERACTION 

The system we want to study is made of two clouds, 
i = 1,2, of dipolar Fermi gas in the normal phase confined 
by the potential Eq. ((T|). Moreover we assume that the 
dipoles are oriented along the z axis by an electric field. 



Nk 



exp 



+ y 2 + k 2 (zt L)"' 



(8) 

where W± and W z — W±/n are the widths of single 
cloud and 2L is the distance between the clouds, which, 
for our parameters, namely a strong confinement in the 
z direction, is very close to the distance 2z$. 

Inserting the gaussian anzatz Eq. (JSJ) into the single 
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TABLE I: The cloud size for N = 22 40 K 87 Rb molecules 
with ujz/2-k — 10 kHz (corresponding to a z = 8.89 x 10 _6 cm) 
and dipole momentum d = 0.56 D. 



X 



W±, 10" 



1.8 
2.9 



W„ 10" 



2.05 
1.54 



TABLE II: The size of cloud for classical gas of 40 K 87 Rb at 
T = Tf with dipole momentum d = 0.56D and trapping 
frequency cj z = 10 kHz (a z = 8.89 x 10" 6 cm)). 



A 


W±, 10- 4 cm 


W z , 10- b cm 


10 


2.8 


2.8 


20 


4.5 


2.25 


40 


7.15 


1.79 




FIG. 3: Frequency shift of ui out /uj±_ for the out-of-phase mode 
for a deg 
Table U 



for a degenerate gas of 40 K 87 Rb, with parameters as given in 



cloud energy one obtains 



E(W x ,k)= ( - 



5/2 



h 2 N f&TT 2 Nn 
2tott V W 3 



mNu 2 Wf ( 1 + 



_*1 

2k 2 



± 

d 2 N 2 K 
Sy^W 3 



2/ .3 



m, (9) 



where we introduce the standard notation (see e.g. [la ]) 
1 + 2k 2 3k 2 arctan y/l — n 2 



/(«) 



1-2k 2 



(l_ K 2)3/2 



Once the densities are known we are in the position 
to calculate the frequency uj ou t- The problem has many 
parameters and to be concrete we consider a gas of bi- 
atomic molecules 40 K 87 Rb and reasonable experimental 
values for the number of molecules and for the trapping 
potentials (see Table U). The result for the out-of-phase 
mode frequency as a function of the clouds' distance and 
for different value of A are reported in Fig. [3J 

We see that for small enough distance the larger the 
cloud (larger A for fixed N) the smaller the effect. This 
can be easily understood in terms of the potential of a 
single disk of radius W± on a probe dipole. Indeed at a 
distance z <§C W± the potential decays with W± , which is 
a general result independent of statistics. On the other 
hand we have that the asymptotic behavior of the fre- 
quency shift at large distance is oc <J1 + C/L 5 with C a 
constant, which is the result one immediately obtaines 
by considering just two trapped dipoles. It can also 



be easily shown that considering spherical clouds with 
W± = W z — W in Eq. ([7]) the frequency of out-of-phase 
dipole mode reads 



w_l 1- 



V2Nd 2 h(L/W)\ 



1/2 



(10) 



where %) = e - 2y2 (4y 5 + 6y 3 + 9/2 2 /)-9/2 1 /f Erf(\/2y), 
which approaches a constant for large values of y. 

The results reported in Fig. [3] show how the frequency 
shift can be large enough, for the chosen parameters, to 
be experimentally measurable. Moreover we checked that 
the results are the essentially the same using Thomas- 
Fermi profiles, instead of the gaussian ansatz Eq. ((SJ. 

It is useful to compare the above predictions for the fre- 
quency shifts calculated for a zero temperature Fermi gas 
with the ones holding for a Bose-Einstein condensed gas 
or for a classical thermal configuration. To this purpose 
one can still use Eq. (JTJ) with the proper density profiles 
(an inverted parabola for a BEC gas and a Boltzmann 
distribution for a classical gas) . As emerges from Eq. J7]) 
the effect is amplified for smaller radial sizes where the 
gradient of the density is larger. One then understands 
that a Bose gas interacting with a moderate value of the 
scattering length will provide larger shifts with respect to 
both a Fermi gas and a thermal configuration. The shifts 
for a Bose gas were investigated in [1 21 ] where, however, 
magnetic dipolar atomic gases were considered, which 
are characterized by a significantly small value of the 
dipolar coupling constant d 2 in Eq. (U). At present the 
more promising perspectives for realizing electric dipolar 
molecules, where the effect of the dipolar interaction is 
particularly strong, concern the fermionic species which 
ensure better stability conditions and which are already 
available in the thermal regime. 

In Fig. 2] we report the predictions for the frequency 
shifts exhibited by a thermal configuration calculated at 
the temperature T = T F , where T F = Huj ± (6NX) 1 / 3 /k B 
is the Fermi energy and ks the Boltzmann's constant. 
We used simply the gaussian density profiles of Eq. © , 
but with the radii given by the Boltzmann expression 
W]_ = 2ksT/ (muj 2 _) and k = A. The corresponding pa- 
rameters are given in Table [TTJ Comparison with the 
predictions reported in Fig. [3] shows that the effect, for 
the same trapping conditions and number of particles, is 
indeed smaller than for a degenerate Fermi gas, since the 
thermal radii are larger and the densities smaller than 
the ones of the degenerate configuration. 
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IV. CONCLUSIONS 

We have proposed a drag experiment between two non- 
overlapping atomic/molecular clouds (see Fig. [T]and Fig. 
[2J to test the long-range nature of the dipolar potential. 
The method is independent of quantum statistics and 
holds for both degenerate and thermal gases. This ef- 
fect corresponds to the trapped version of the famous 
Coulomb drag exhibited by electrons in uniform bilayer 
systems. The realization of such a drag experiment would 
provide a direct and easy signature of the long-range na- 
ture of the dipole interaction. 



FIG. 4: Frequency shift of Ld ou t/uj± for the out-of-phase mode 
for a classical gas of 40 K 87 Rb at temperature T = Tf, where 
Tf is the Fermi temperature of the gas. Parameters are given 
in Table HQ. 
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